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CONTINUITY AND STRICT POSITIVITY OF THE MULTI-LAYER 
EXTENSION OF THE STOCHASTIC HEAT EQUATION 

CHIN HANG LUN AND JON WARREN 


Abstract. We prove the continuity and strict positivity of the multi-layer extension 
to the stochastic heat equation introduced in [OWll] which form a hierarchy of 
partition functions for the continuum directed random polymer. This shows that 
the corresponding free energy (logarithm of the partition function) is well defined. 
This is also a step towards proving the conjecture stated at the end of the above 
paper that an array of such partition functions has the Markov property. 


1. Introduction 

In [OWll] O’Connell and Warren introduced the following: for each n = 1,2,..., 
t > 0 and X, y gR define 

Znit,x,y) =pt{x-y)^(l + '^ [ [ Rk{s,y';t,x,y)W®’'{ds,dy')\ ( 1 ) 

where Afe(t) = {0 < si < S2 < • • • < s* < t}, s = y' = {y\,...,y'^ and 

Rk{s,y';t,x,y) is the fc-point correlation function for a collection of n non-intersecting 
Brownian bridges each of which starts at x at time 0 and ends at y at time t. pt{x — y) 
is the heat kernel The integral is a fc-fold stochastic integral with 

respect to space-time white noise, see Section 2 for the definition of such integrals. It 
was shown in [OWll] by considering local times of non-intersecting Brownian bridges 
that the infinite sum in the definition is convergent in with respect to the white noise. 

Observe that u = Zi is the solution to the (multiplicative) stochastic heat equation 
(SHE) with delta initial data: 


dtu{t, X, y) = (| Ay -f W{t, i/)) u{t, x,y), t e {0,oo),y e R, 
m(0, X, y) = 5{x — y), x G R. 


By a solution to the above we mean a random field u which satisfies almost surely the 
mild form 


u{i,x,y) =pt{x-y) + 


Pt-siy 


y')u{s,x,y') lU(ds,dy')- 


( 3 ) 


Jo Jr 

Iterating equation (3) multiple times gives the chaos expansion (1) for n = 1. One can 
express the solution u{t, x, y) in a more suggestive notation: 


u{t,x,y) = pt{x 


y)K,y,t 



W{s, bs) ds 


( 4 ) 


where 6 is a Brownian bridge that starts at x at time 0 and ends at y at time t and 
denotes the corresponding expectation, ^xp is the Wick exponential defined by 


(g’xp(Mt) :=exp (Mt - ^{M, M)t), 
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for a martingale M. The Feynman-Kac formula (4) is not rigorous as it is unclear how 
one would define the integral of the white noise along a Brownian path and moreover 
to exponentiate such an expression. However, Taylor expanding the exponential, then 
switching the expectation with the infinite sum and evaluating the expectation, one 
obtains the chaos expansion of u. With this in mind, (4) can be thought of as a short 
hand for the chaos expansion (1) in the case n = 1. On the other hand, one can obtain 
an rigorous expression by replacing W in (4) with a smoothed version of the space-time 
white noise. Indeed, Bertini and Cancrini showed in [BC95] that such expression has a 
meaningful limit as one takes away the smoothing and that the limit solves the SHE. 
With this Feynman-Kac interpretation, one can think of the solution to the stochastic 
heat equation as the partition function (up to a multiplication by the heat kernel) of 
the continuum directed random polymer [AKQ14a]. 

Analogously, we write 


Zn{t,x,y) 


nmX 




'x,y,t 


^xp 


E 


W{s,Xl) 



( 5 ) 


where ..., W", 0 < s < t) denotes the trajectories of the above mentioned collection 
of n non-intersecting Brownian bridges and corresponding expectation. In 

the same manner as in the n = 1 case, (5) should be thought of as the short hand for 
the chaos expansion (1). Therefore, in view of (5) one can interpret as the partition 
function (up to a factor of the heat kernel) of a natural extension of the continuum 
directed random polymer involving multiple non-intersecting Brownian paths. 

Since the work of Bertini and Giacomin [BG97], it is widely accepted that the loga¬ 
rithm of u is the Cole-Hopf solution to the KPZ equation [KPZ86], 

dth{t, x) = d‘lh{t, x) + {dxh{t, x)Y + W{t, x), (6) 

with narrow wedge initial condition. This solution arises as the scaling limit of the 
corner growth model under weak asymmetry. The Gole-Hopf solution to the KPZ equa¬ 
tion via the Feynman-Kac formula (4) can be seen as the free energy of the continuum 
directed random polymer. With this interpretation the Gole-Hopf solution can be re¬ 
garded as the continuum analogue of the longest increasing subsequence of a random 
permutation, length of the first row of a random Young diagram, directed last passage 
percolation and free energy of a discrete/semi-discrete polymer in random media etc., 
see [BDJ99a], [BDJ99b], [BOOOO], [Joh99], [JohOla], [PS02], [JohOS], [GOSZ14] and the 
references therein. In each of these discrete models, there is further structure provided 
either by multiple non-intersecting up-right paths on lattices, multi-layer growth dy¬ 
namics or Young diagrams constructed from the RSK correspondence. The work in the 
above mentioned references have shown that in some cases, utilisation of this additional 
structure have lead to derivations of exact formulae for the distribution of quantities 
of interest. The above mentioned discrete models provide examples of what is called 
integrability or exact solvability. The motivation for introducing the partition functions 
Zn, which are the continuum analogue of the structures mentioned above, is that they 
should provide insight to the integrable structure in the continuum setting. 

The main result of this paper is that the continuum partition functions possess some 
nice regularity properties. 


Theorem 1.1. For all n > 1, the function {t,x,y) i—>■ Zn{t,x,y) has a version that is 
continuous over (0, oo) x R x K. Moreover, 

P[E„(t, x,y) >0 for all t > 0 and x, y G R.] = 1. 
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Now define for n = 1, 2,... 


hnit,x) = log 



( 7 ) 


with the convention that Zq = 1, then hi{t,x) is the Cole-Hopf solution to the KPZ 
equation with narrow wedge initial data. An immediate corollary to the above theorem 
is 

Corollary 1.2. For all n > 1, is well defined and it is a continuous function of (t,x) 
over (0, oo) x R. 

The collection {/i„, n > 1} represents a multi-layer extension to the free energy of the 
continuum directed random polymer. It is the analogue in the setting of the KPZ of the 
multi-layer PNG or its discrete counterpart studied in [PS02] and [Joh03] respectively. 

We mention here the work of [CH13]. The authors showed the existence of a collection 
of random continuous curves such that the lowest indexed curve is distributed as the 
time t Cole-Hopf solution to the KPZ with narrow wedge initial data. It is believed (see 
[CH13, Conjecture 2.17]) that for each t > 0 hxed, their collection of curves is equal to 
{hn{t,x) : n > l,x G K} defined by (7). Proving this will give an alternative proof of 
the continuity and strict positivity of Zn at a fixed time t. In this paper, we provide a 
direct proof of this and furthermore our proof gives a stronger result since t can vary 
over (0, oo). 

The continuity and strict positivity oiu = Zi was proved by considering its mild form 
which suggests that to prove Theorem 1.1 one could consider the evolution equation 
satisfied by Z„. By considering a smooth space-time potential, the authors in [OWll] 
showed that Zn should satisfy a certain SPDE, see [OWll, Proposition 3.3 and 3.7], 
however unfortunately it is not immediately obvious that this SPDE makes sense in the 
white noise setting. Instead, we shall show that a natural extension of does satisfy a 
rigorous evolution equation which can be regarded as a multi-dimensional stochastic heat 
equation. This allows us to derive the continuity and strict positivity of the extension 
and from which Theorem 1.1 follows as a corollary. 

Denote by Wn the Weyl chamber {x G R" : xi > X 2 • • • > x„}, then for n = 1, 2,..., 
t > 0 and X, y G Wn define 



where Rk is the fc-point correlation function of a collection of n non-intersection Brow¬ 
nian bridges which starts at x at time 0 and ends at y at time t. p*(t,x, y) = 
det[pt(xi — is by the Karlin-McGregor formula [KM59] the transition density 

of Brownian motion killed at the boundary of Wn- It was proved in [OWll, Proposition 
3.2] that Kn also satisfies a Karlin-McGregor type formula: 


Kn{t,x,y) = det[u(t,x*,2/j)]”j^i. 


( 9 ) 


where each term in the determinant are solutions to (2) each driven by the same white 
noise. Now, define for t > 0, x, y G W^ 



( 10 ) 


where A(x) = ni<i<j<n(®i ~ Vandermonde determinant. It follows from (8) 

that Mn has chaos expansion 



1 + X! / f 7?fe(s,y';t,x,y) IE®'"(ds,dy') 


JAkit) JR’- 


( 11 ) 
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By (9) and the continuity of the solution to the stochastic heat equation, it is easy 
to see that y) is almost surely continuous on (0,t) x Wn x Wn and is zero 

on the boundary of Wn x Wn- It follows that y) is continuous in the interior 

^ '^n- By [BBO09, Lemma 5.11], p*(t,x,y)/A(x)A(y) is a smooth function of 
(x, y) over R." x R" and since the fc-point correlation function Rk extends continuously 
to the boundary of the Weyl chamber, see Section 2.4, we see from its chaos expansion 
(11) that M„(t,x, y) is defined for x, y S dWn- This also suggests that M„(t,x, y) is a 
continuous function on Wn x Wn- Furthermore, from (9) we see that M„ being a ratio 
of determinants is a permutation symmetric function of its spatial variables, that is for 
any permutations tt, ct of {1,..., n}, M„(t, ttx, ay) = Mn{t, x, y). Hence, we can extend 
Mn by symmetry to a function on R" x R" and we will show that there exists a version 
of Mn that is almost surely strictly positive and continuous on the whole of R" x R" 
and for all t > 0. Moreover, when all the x coordinates are equal and likewise for y, 
agrees up to a multiplicative constant with Z„, that is 


Mn{t,al,bl) = Cn,tZn{.t,a,b), (12) 

where Cn,t ■= and 1 = (!,...,!). Equation (12) was shown 

to hold in [OWll] but there the continuity of Mn on the boundary of Wn was only 
established in an sense; here we extend it to almost sure continuity. Note that (9) 
suggests that Er„(t,x, y) and M„(t,x, y) can be regarded as the stochastic analogue 
of p*(t,x, y) and p* (t, x, y)/A(x)A(y) respectively where the latter has limit at the 
boundary equal to Cn,tPt(a — &)"• 

In Section 4, we will show that for all (t, x, y) S (0, oo) x R” x R", M„(f, x, y) satisfies 
almost surely the mild equation 


M„(f,x,y) ^ +An[ [ Qt_s(y,y')^n(s,x,y') dy', IT(ds,dp'i) 

=: J„(t,x,y)+/n(t,x,y), (13) 

where An = l/(n — 1)! is a combinatorial constant, dy(. = dp 2 ■ • ■ dy„ and 


Qt(x,y) = -^^p*n{t,x,y), 

is the transition density of Dyson’s Brownian motion starting from x G Wn and ending 
at y £ Wn and it satisfies 


Qt{al,y) = c„,tA(y)^ Y[pt{yi - a)- (14) 

i=l 

We can extend Qt by symmetry to a function on R” x R" and so the integral over R” 
in the mild equation (13) is defined. 

Consider also the following integral equation for (t,y) £ (0,oo) x R", 


Mn{t,y) 


= ^ [ 5(y')0i(y:y') dy' 

n' J-g^n 

+ An [ f Qt-s{y,y')Mnis,y') dy'^W{ds,dy[) 
Jo Jr^ 

=: Jn{t,y) + In{t,y), 


(15) 


where g : R" —> R is permutation symmetric and may be random but independent of 
the white noise. The function g is the initial condition for equation (15) in the sense 
that 


lim — 

(->■0 n! 


9{y')Qt{y,y') dy' = lim f g{y')Qtiy,y') dy' = p(y). 
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On the other hand, we say that y) is the solution started from a delta initial 

data at x even though strictly speaking it is the ratio of itr„(t,x, y), which can be 
shown to satisfy an integral equation similar to (15) with delta initial condition, and 
the product of Vandermonde determinants A(x)A(y). To emphasise the initial data we 
sometimes write M®(t,y) instead of M„(<,y). 

We now state the main results regarding the solutions of (13) and (15) from which 
Theorem 1.1 follows as a corollary by (12). Let be the collection of Borel measur¬ 

able subsets of R with finite Lebesgue measure and let W = (Wt{A),t > 0, A G .^f,(R)) 
be space-time white noise on a complete probability space (fl,^, P) endowed with a 
right-continuous filtration {^t)t>o such that W is ^(-adapted and Wt{A) — Ws{A) is 
independent of for all A G ^bO^)- From now on we fix this filtered probability space 
(11, P)- We use E to denote the expectation with repect to P and for p> 1, 

II • lip = (E[| • |P])^/P denotes the L^’(H) norm. Throughout this paper, Cp < 2y/p is the 
constant appearing in the Burkholder-Davis-Gundy inequality. 

Theorem 1.3. (a) Suppose that g is ^Q-measurable and symmetric and satisfies 

for all p > 2, supygRn ||ff(y)||p < Kp,g < oo> then there exists a solution 
[Mn{t, y), (t, y) G [0, oo) xM") to the integral equation (15) that is unique (in the 
sense of versions) in the class of all random fields (u(t, y), (t, y) G [0, oo) x M") 
that satisfy sup(( y)g[g lk(^!y)llp < oo for all T > 0. The solution satisfies 

for all p>2 



(16) 


for a constant A > 0 depending on n. 

Moreover, M„ has a version such that (t,y) i—>■ M„(t,y) is locally Holder 
continuous on (0, oo) x M" with indices a <1/2 in space and a < 1/4 in time. 
(b) There exists a unique solution (M„(t,x, y) G (0, oo) x M" x R") given by the 
chaos expansion (11) to the integral equation (13) such that for all p > 2 and 


t > 0 


sup ||M„(t,X,y)||2 < C^^pt 


(17) 


x.yGR” 


for some constant Cn,p > 0. 

Moreover, Mn has a version such that (t,x,y) i—>■ Mn{t,x,y) is locally Holder 
continuous on (0,oo) x R" x R" with indices a < 1/2 in space and a < 1/4 in 
time. 

Theorem 1.4. Let g be as in Theorem 1.3(a) with the additional property that g is 
non-negative almost surely and P[(jf(y) > 0 for some y G R"] = 1. Then the solution 
to (15) satisfies 

P[M®(t, y) > 0 for all t > 0 and y G R"] = 1. 

Let Mn be the random field defined by (11) then 

P[M„(t, X, y) >0 for all t > 0 and x, y G R"] = 1. 

Comparing (13) and (15) with (3), we see that they have a similar form to the mild 
equation of the SHE which has been well studied. It has been shown for various initial 
data that the solution is Holder continuous with indices up to 1/2 in space and up to 
1/4 in time. For example, the case with a bounded initial data was studied by Walsh in 
[Wal86] . Bertini and Cancrini stated the Holder continuity in [BC95] for a class of initial 
data which includes a delta function. More recently, Chen and Dalang [CD14a] proved 
the Holder continuity for a non-linear SHE with initial data p being a signed Borel 
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measure over K. such that (|^| *pt){x) < oo for all t > 0 and a; G K.. For other variants 
of the SHE see for example [CJKS14], [Shi94], [SSS02] and the references therein. 

In each case the tool used to prove the continuity of the solution is Kolmogorov’s 
continuity criterion. Denote the stochastic integral term of (3) by I(t,y) then the key 
is to show that 

E[| J(t, y) - I{t', yT] < C{\y - y’Y'^ + \t - t'\P/^), 

for p large enough. This in turn requires showing some continuity estimate for the 
heat kernel and in our case, estimates for the kernel Qt, see Theorem 3.1 below. These 
estimates get increasingly involved for increasingly less regular initial data due to the 
pth moments E[|u(t,?/)| p] of the solution being unbounded as t I 0 or as y —>■ oo or 
both. However for certain initial data such as a delta function, even though the pth 
moments blow up as time t 0, they are for any fixed positive times uniformly bounded 
in space and thus one can in effect isolate the effects of the initial data by solving the 
equation for a small time and then start afresh with the current solution as the new 
initial condition. This is the case with M„(t,x, y). We will show that for all positive 
times t, E[|M„(t,X, y)|P] is bounded uniformly in space for all p which puts us in the 
situation of (15) with g having uniformly bounded pth moments for which continuity is 
easier to obtain. 

The strict positivity of the solution to the stochastic heat equation was first proved 
by Mueller in [Mue91]. He showed that if the initial data / is non-negative, continuous 
with compact support with f{x) >0 for some a: G K., then for alH > 0 

P[u(t, x) > 0 for every x G R] = 1. 

Bertini and Cancrini proved a weak comparison principle using the Feynman-Kac for¬ 
mula and used it to extend Mueller’s result to a delta type initial data. Shiga in [Shi94] 
proved the stronger statement 

f‘[u{t,x) > 0 for every x G K and every t > 0] = 1, 

for initial data being continuous function such that the tails grow no faster than 
for all A > 0. More recently, Moreno Flores in [Flol4] proved the strict positivity of 
the solution for delta initial conditions, using a convergence result of a discrete poly¬ 
mer model to the SHE, see [AKQ14b]. Chen and Kim [CK14] further generalised the 
strict positivity result to the fractional SHE, which includes as a special case the SHE 
considered here, for measure-valued initial data by adapting Shiga’s method. 

In all of the proofs above (except for the polymer proof) a key result is a large 
deviation estimate on the stochastic integral term of the solution. Mueller proved such 
result using the fact that integrals of the type f{s, y) IF(ds, dy) can be considered 
as a time-changed Brownian motion. Chen and Kim using a method of [CJK12] derived 
a similar estimate for the fractional SHE using Kolmogorov’s continuity criterion. We 
will adapt the approach of [CK14] since we will first derive the necessary estimates in 
order to prove Holder continuity anyway. 

The outline of the paper is as follows. In Section 2.1 we first briefly recall integration 
with respect to space-time white noise and multiple stochastic integrals. In Section 2.2 
we derive an upper bound on the L^(VL) norm of stochastic integrals which will be used 
repeatedly in this paper and we discuss briefly non-intersecting Brownian bridges in 
Section 2.4. We then prove some estimates on the transition density Qt in Section 3 
which are central to the proof of existence and continuity. The existence, uniqueness 
and moment estimates part of Theorem 1.3 will be proved in Section 4. The proof of 
Holder continuity is in Section 5. Finally, in Section 6 we prove a strong comparison 
principle for the integral equation (15) of which Theorem 1.4 is a corollary. 
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2. Preliminaries 

2.1. White Noise and Stochastic Integration. In this section we briefly recall the 
Walsh stochastic integral with respect to white noise, see for example [Wal86], [Kho09] 
and [Dal99] for details. Let be the collection of Borel measurable subsets of 

with finite Lebesgue measure. A white noise on is a mean zero Gaussian random 
field with covariance function 

E[W{A)W{B)] = \An B\, for all A,B G 

where | • | denotes the Lebesgue measure on R.'^. We will only consider the case d = 2 and 
we interpret one of the dimensions as time. More precisely, we define a space-time white 
noise {Wt{A)A > 0,A G by Wt{A) := lT([0,t] x A) on a filtered probability 

space (H, as described above Theorem 1.3. 

A random field / is elementary if it is of the form 

f{s,y) = Al(a_b](s)lA( 2 /), 

where X is bounded and .^a-measurable and A G .^(R). A simple function is a finite 
linear combination of elementary functions. We say that a random field / is predictable 
if it is measurable with respect to the cr-algebra generated by the simple functions 
and we say that / G 3^2 if it is predictable and / G x [0,oo) x R). According to 
Walsh’s theory, [Wal86], {Wt{A)} belongs to a suitable class of integrators called worthy 
martingale measures and the integral 

f [ fis,y)W{ds,dy), 

Jo Jk. 

is defined for all / G ^^ 2 - 

Now we turn our attention to multiple stochastic integrals which appear in the chaos 
series in the introduction. Let fc > 1 and let / G L|([0, x R'^) such that /(tts, Try) = 
/(s,y) for all (s,y) G [0,t]^ x R^ and tt G Sk where Sk is the set of permutations of 
{1,..., fc} and TTS = (s^i,..., s,rfe)- Let Ai,... ,Ak be disjoint subsets of [0, t] x R. An 
elementary function in L|([0,t]^ x R^) is a function of the form 

k 

/(s,y) = n n ; y-TTi) € (1^) 

For such / we define the fc-fold integral by 

P k 

/ /(s,y) W®'=(ds,dy) = fc!TTw(A,). 

1,4]'= dR'= 

It can be shown that linear combinations of functions of the form (18) are dense in 
L|([0,t]^ X R^) and that for an elementary /, the integral (/ • W)k satisfies an ltd 
isometry, hence for a general / G L|([0,t]^ xR^), we define {f ■W)k = lim„_>.oo(/n-lE)fc 
where {fn}n>i is a sequence of elementary functions such that /„—>■/ in L^([0, t]'' xR^). 
The resulting integral is a mean zero random variable with covariance given by 

E[(/ ■ W)k{t){g ■ W)k{t)] = (/,5)l2([o.4]'=xr'=)- 


if-W)kit)= [ 


(19) 
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For / G tY X that are not symmetric, we define its integral by first symmetris¬ 

ing / via 

/(s,y) '■= 

TrGSfc 


and then define 

if-wut) = if-wut). 

Finally, for functions / defined on Afc(t) x R.*^, for example the /c-point correlation 
function Rk appearing in (1) and (8), we first extend it to a function on [0, by setting 
it to be zero for s ^ Afe(t) and then define 


f [ /(s,y) W^®'=(ds,dy) 
JAk{t) JW’ 


if-wut). 


Now define a time reversed white noise W by fF([0, s] x A) = W{[t — s, t] x A), s <t 
and A G S§b{^)- We will need the following result for the proof of continuity in the 
initial data. 


Lemma 2.1. Let f G L|([0,t]*^ x R^) then 

f [ /(s,y) W®'=(ds,dy)= / [ /(t-s,y) W®'=(ds,dy) a.s.. 


no,t]’’ 

where t — s = {t — si,... ,t — Sk)- 


Proof. The result in the case when / is an elementary function of the form (18) follows 
from the definition of the integral and the definition of W. For general / G L|([0, x 
R^), let {fn}n>i be a sequence of elementary functions converging to /. The result of 
the lemma holds for (/„ • W)k{f) for all n and by taking limits we see that the result 
also holds for (/ • W)k{t). □ 


2.2. PP Bounds on Stochastic Integrals. The following estimate is a useful bound 
on the LP{n) norm of stochastic integrals; it can be considered as a version of [CK12, 
Lemma 2.2] or [FK09, Lemma 3.3] adapted to the present setting. Recall that for brevity 
we denote y(. = di /2 • • ■ dy], and Cp < 2y/p is the constant appearing in the Burkhoider- 
Davis-Gundy inequality. 

Lemma 2.2. Define a random field {f{t,y);{t,y) G (0, 00 ) x R") by 

f{t,y)=[ [ rt_^(y,y')w(s,y') dyl kF(ds,d?/'i), 

Jo JR" 

for a suitable random field w and Tt(jj,y') is a non-random measurable function on 
(0,oo) X R" X R" such that Ft_s(y, y')w(s, y') dy(. is integrable in the sense of 

Walsh for all {t, y) G (0, 00 ) x R". Then for all integers p > 2, t > 0 and y G R" 

||/(Ly)llp < Cp / [ ( [ rt_,(y,y')||w(s,y')||p dy',^ dy'^ds. 

Jo Jr V./r"-i / 


Proof. Fix t and y, then by the Burkholder-Davis-Gundy inequality applied to the 
martingale (/J’/r„ Tt-s{y,y')w{s,y') dy(. W{ds,dy[), r G [0,t]), we have 


\\f{fy)\\l<cl 


rt-s(y,y')^c(s,y') dy(, 


2 

dy'ids 


p/2 
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Applying Minkowski’s integral inequality [Kal02, Corollary 1.30] twice, we obtain 


\\f{t,y)\\l<cl f [ f rt_s(y,y')w(s,y') dy' 
Jo Jr Jr"-i 


dj/]ds 


<c; 


f f ( f rt_s(y,y')||'!T(s,y')llp dy',^ dy'^ds, 
Jo Jr \Jr"-i / 


as required. 

Lemma 2.3. For all k > 1 and f € L^(Ak(t) x R^) we have 


□ 


'Afc(t) Jr>^ 


/(s,y)M^®'=(ds,dy) < 


[ [ /(s,y)^dyds. 

JAkit) Jr'‘ 


Proof. Since multiple stochastic integrals on Ak{t) coincides with iterated stochastic 
integrals, applying Burkholder-Davis-Gundy inequality and Minkowski’s integral in¬ 
equality k times gives the desired upper bound. □ 

2.3. Predictability of Random Fields. Recall that the Walsh integral is defined for 
random fields in ^ 2 , see Section 2.1 above, therefore it is convenient to have a set of 
conditions to verify the predictability of a random field. The following result is from 
[CD14b, Proposition 3.1] which is an extension of [DF98, Proposition 2] to space-time 
white noise. 

Proposition 2.4. Let t > 0 and suppose a random field {f{s,y),{s,y) € (0,t) x R.) 
satisfies 

(i) / is adapted, that is for all {s,y) G (0,t) x K, f{s,y) is ^s-measurable; 

(ii) for all (s,?/) G (0,t) x R, ||/(s,2/)||2 < 00 and {s,y) f{s,y) is 
continuous on (0,t) x R; 

iii) IoIr \\fis,y)\\j dyds < 00 . 

Then /(•, •)l(o.t)(’) ^ -^2 and 

f [ fis,y)W{ds,dy), 

Jo Jr 

is a well-defined Walsh integral. 

In the sequel we will need to integrate functions of the form: for some random field 
M, let f{s,y[) = Qt-s{y,y')M{s,y') dy*. (Note that we have suppressed the 

dependency of / on t and y to keep the notation simple). The following proposition 
provides convenient conditions to verify the integrability of such a random field. 

Proposition 2.5. Let t > 0 and y G R". Suppose the random field (M(s,y'), (s,y') G 
(0, t) X R") satisfies 

(i) M is adapted i.e., for all {s,y') G (0,t) x R", M{s,y') is -'measurable; 

(ii) {t,y') H> M(s,y') is L'^(LI)-continuous on (0,t) x R”; 

(iii) sup(^_y,)g(o_t)xR„ |lM(s,y ')||2 < 00 ; 

Then (f(s,z),(s,z)€(0,t)xR) defined by f{s,y[) = Qt-s{y,y')M{s,y') dy',, is 

in J ^2 and 

r f f{s,y[)W{ds,dy[), 

Jo Jr 

is a well-defined Walsh integral. 
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Proof. We will show that / satisfies the three assumptions of Proposition 2.4. Since 
Qt-s{y,y') is continuous and deterministic, Qt-s{y,y')M{s,y') is adapted by (i) and 
so the integral Qt-s{y, y')M{s, y') dy^ is also adapted. Assumption (iii) of Propo¬ 
sition 2.4 follows from (iii) above since by Lemma 2.2 and Lemma 3.8 below, we have 
for some constant C 



\\fis,y'i)\\l dy'ids 

< 


Qt_^(y,y')||M(s,y')||2 dy(. dy'^ds 


< 


sup \\M{s,y')\\l [ [ ( f Qt_s(y,y') dy',^ 

')6(0,t)xR" Jo JR VJR"-i / 


(s,y 

<20^/"^ sup ||M(s,y') 


2 

dy'ids 


It remains to show the L^(n)-continuity of /. We wish to show that for each (s,y) G 
(0,t)xR, \im.(^u,z)^{s,y)\\f{u,z)-f{s,y )\\2 = 0 . Let /i > 0 and suppose zi £ [y[-h,y[+h] 
and u £ [s/2, {t + s)/2]. Then by the Harish-Chandra formula (22) and equation (23) 
below, we have 


Qt-uiy, z) < Cn{t - n)-"'/2A(z)2 H e-(J.--dV2d-«) 

i=l 


< 


2"V2c 


(t — s)"^/2 




+ \Zi 


( V-i — 2,) 

)2g- 2(t-^/2) g- 


yi-^\yi+h\yi 

2{t-s/2) 


2<i<j<n 


The last line is integrable with respect to dz* = dz 2 • • • dzn and so by the dominated 
convergence theorem, the continuity of Qt and assumption (ii), the right hand side of 

||/(m,2;i)-/(s,2/'i)||2 


1/2 


< sup ||M(s,y )||2 / |Qt_„(y, (zi,z,)) - Qt_s(y, (i/;,z,))| dz, 

(s,y) \ 

+ [ (5t-„(y, (i/i,z*))||M(s, (zi,z*)) - M(s, (y'i,z ,))||2 dz* 


converges to zero as {u, zi) —>■ {s,y[). Finally, an application of Proposition 2.4 com¬ 
pletes the proof. □ 


2.4. Non-intersecting Brownian Motions. Dyson Brownian motion introduced in 
[Dys62] can be realised as the eigenvalues of Hermitian Brownian motion, an n x n Her- 
mitian matrix whose entries are (up to the Hermitian condition) independent standard 
complex Brownian motions. The eigenvalues of such a matrix is a Markov process with 
state space Wn with transition density Qt{x,y). It also arises as the Doob /i-transform 
of Brownian motion killed at the boundary dWn with h{x) = A(x) (see for example 
[Gra99] and [KT07]). 

One can construct bridges of Dyson Brownian motion, which we will call Dyson Brow¬ 
nian bridge or non-intersecting Brownian bridges, using the framework of [FPY93]. For 
X, y £ Wm a collection of non-intersecting Brownian bridges Xt = {X },..., X'f) starting 
at X at time 0 and ending at y in time t is a process whose law is absolutely continuous 
to that of Dyson Brownian motion started at x with Radon-Nikodym derivative equal 
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to 

Qt—s{^sj y) 

Qt(x,y) 

In particular, for 0 < si < ... < Sfe < t, the law of , • ■ •, is given by the density 

Qsi (x, y^) (y'~^ y')Qt-s^, (y^, y) 

Qt(x,y) 


The above is well defined at the boundary of the Weyl chamber by (14); in particular, 
taking limits as x —>■ al, y —>■ &1 where 1 = (1,..., 1) one obtains 

^ A(y^)A(y^) UU y^-\ y^) UU - Vj) 


where = 0"=/ The fc-point correlation Rk appearing in (8) is defined as the sum 
over ii,..., ffe for 1 < V < n, 1 < r < A: of the densities of the process (A*),..., A|^): 





x,yi) 11^2^(51 - g^-l.y' ^ 

K(^x,y) 


r)p*jt 


sk,y'",y) 


k n 


nn-isi 

^=1 17^*1 


Notice that the integrand above is symmetric in the permutation of its arguments 
{y[,... ,yl^) for all 1 < I < k and so we can rewrite each integral over Wn-i as integrals 
over multiplied by a factor of 1/n!. Moreover, by symmetry each term in the 

sum over ii,... ,ik gives the same contribution. There are in total of such fc-tuples 
and hence we can rewrite the correlation function ((si, j/^), ..., (sfe, yf ); t, x, y) := 
i?fc(s,yi;t,x,y), yi = {yl, ... ,yf) as 


A 


k 

n 



K(suX,y^)nL 2 Pn('S» - S^-l,y^ \y")Pn{t 

K(i>x,y) 


sfe,y'',y) 


k 71 

(20) 


i=l i=2 


where A„ := l/(n — 1)!. For each k we have chosen to leave the first coordinate of y^ 
and integrated out the rest but this choice is arbitrary by symmetry. Note that this is 
also the reason for the form of the stochastic integral term in (13). 

In the sequel we will need to bound integrals of the square of the fc-point correlation 
function R^. Correlation functions of densities given by a product of determinants have 
been studied extensively in the context of determinantal point processes, see for example 
[Joh06] and [Borll]. They can be expressed as a determinant of a matrix whose entries 
are given by some kernel function. However for general start and end points x and y 
this kernel function is difficult to compute, but since all we need is the integral of the 
square of Rk it is not necessary to compute Rk explicitly and so we will not pursue this. 
Instead, the next two results proved in [OWll] which expresses the integral of R\ in 
terms of intersection local times of Brownian bridges will be used. Let X = (A^,..., A") 
and Y = (F^,... ,F") be two independent copies of a collection of n non-intersecting 
Brownian bridges which start at x at time 0 and end at y at time t and let denote 

the corresponding expectation of the joint law of the bridges. Let Lt{X'’ — Y^) be the 
local time at 0 of the difference A* — X^. Then we have 


Lemma 2.6. Fix n > 1. For all integers k>\ and all t > 0, x, y G Wn the following 
holds 




Rk{s,y';t,x,yf 


dy'ds = 


1 h^a.y 


71 

[(E u(-v 




k-\ 
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The following is used to bound the above moments of local times. 


Lemma 2.7. For all a > 1 and t > 0, there exists a constant C > 0 such that 


The above two lemmata shows that for each t > 0, \\Zn(t,x,y )\\2 < oo uniformly in 
X and y and thus the chaos series (1) is convergent in The same is also true for 

(8). We point out here that the bound on the pth moments of Mn{t,x,y) can in fact 
be written as 


||M„(t,x,y)||2 < 2 


VA(x)A(y) 


E: 


X,Y 

x,y;i 


exp (2cl 

i,j=i 


( 21 ) 


The bound (17) in Theorem 1.3(b) then follows from the above by Lemma 2.7. 


3. Estimates on Qt 


From now on we drop the bold typeface for vectors in R" or Wn since we will only be 
working with functions of multi-dimensional spatial variables so there is no longer any 
risk of confusion. 

Before proving Theorem 1.3 we need estimates on various quantities involving the 
kernel Qt- The following known as the Harish-Chandra/Itzykson-Zuber formula [IZ80] 
provides a useful alternate expression for Qt'. 

= c„ / exp (Tr YUXU^) dU, (22) 

A{x)A{y) Ju[n) 

for Hermitian matrices X and Y with eigenvalues Xi,... ^Xn and yi, ■ ■ ■ ,yn respectively. 
Cn = (nr=i^*0 integral is with respect to the normalised Haar measure on 

the unitary group lA{n). Furthermore, the integrand above is bounded uniformly in U 
as the following bound from [MRTZ06, Lemma 1] shows 

sup exp ( --Tr(Zly - 17 L>xC7'I')^') < (23) 

UGW(ti) ^2 ^ tXl 

As mentioned in the introduction, Qt(x, y) is well defined on the boundary of the Weyl 
chamber and since it is a product and ratio of determinants, it is permutation invariant 
and so we can extend Qt to a function on K" x K" by symmetry. Denote Kt{x, yi) := 
/gn-i Qt(a^, y) nr=2‘^2/* X := Ki. The following result strongly indicates the 
continuity of M„; in fact it is a key estimate in its proof in Section 5. 


Theorem 3.1. (a) There is a constant Ci > 0 depending only on n such that for 

all t > 0 and x, z G R." we have 

f f {Ks{x,y) - Ks{z,y)y dyds < Ci\x - z\, 

Jo Js. 

(b) for all a < 1/2 and T > 0 there are positive constants C 2 '■= 02 ( 0 :,n,T) and 
C'a := Coin) such that for all t, u with 0 < u < t < T and x G R", we have 

f f {Kt-u+six,y) - Ks{x,y)Y dyds < C2\t - u\°‘, 

Jo Jr 



Kt-s{x,yf dyds < Colt - 


and 
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The theorem is a consequence of the series of results below. First observe that Qt 
has the following scaling property: 


= t '^^'^Qi{xlVi,y/Vi). (24) 


^ A{x/Vi) 

The left hand side of the inequality in Theorem 3.1(a) is bounded above by 

pOO r / P \ 2 

/ / ( / Qs{x,y)-Qs{z,y)\\yi] dyids 

Jo Jr \ JR"-i / 

^ lo ^ ) - Qi(2:/\/s,y0 dy'ids, (25) 


where we have changed the integration region to [0, oo) in the time integral which results 
in an upper bound due to the positivity of the integrand. The equality follows from the 
scaling property (24) and a change of variables. Theorem 3.1(a) follows from (25) and 
Lemma 3.2 below. 


Lemma 3.2. Suppose a function R{x, y) : M" x R —>■ R satisfies for some constants ci, 

C2 > 0 

J {R{x,y) - Riz,y))‘^ dy <inm{ci,C 2 \x - (26) 

for any x, z € R”, then 

POC 1 P 

/ {R{x/Vi,y)-R{z/Vi,y))‘^ dydt<C\x-z\, 

Jo yt Jw 

with C = dy/cicf- 
Proof. 

[ [ {R{x/Vi,y) - R{z/Vi,y)y dydt 

Jo yt Jr 

P^\x-z\‘^ POO 

< -^dt+ -^\x - z\'^ dt = C\x - z\. 

Jo Vt J^\x-z\^ R' 

□ 


Thus, we need to show that Kix^ y) satisfies the hypothesis of Lemma 3.2. Using the 
inequality (a + 6)^ < 2(a^ + 6^), the left hand side of (26) with K in place of i?, can be 
bounded by 

sf / K{x,yf dy+ [ K{z,yf dy^ < 4 sup \\K{x,-)\\l-,(^^ y 
\Jr Jr J a:GR" 

On the other hand, let r(p) : [0,1] —>■ R", r(p) = (1 — p)x + pz be a parameterisation of 
the straight line from a; to z, then 

K{x,y) - K{z,y) = [ WK{r{p),y) ■ r'{p) dp, 

Jo 

where the gradient is with respect to the first variable of K{-, •) and u ■ v denotes the 
usual inner product of two vectors in R". Then by Minkowski’s integral inequality and 
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Cauchy-Schwarz inequality we have 

- K{z,y)Y dy^ </ II •) • r'(p)||L2(d,/) dp 

< sup |||V-ft:(r(p),-)|||i2(d„) | 2 ;-^l- 

pG[0,l] ^ 

Therefore, in order to verify the hypothesis of Lemma 3.2 we need to show that 

sup / K{x,y)'^ dy < oo, 

kGR" Jm 


and 


sup / |Vi4r(a::, i/)p dy < oo. 

xgK" Jr 

We first concentrate on (28). It suffices to show that 

f 9K , , 2 j 

dy < 00, 

kgR" Jr 


for all j = 1,..., n. Clearly, 


f dK, ,2 , fdK, \ f 

/ dy<sup —(x,y) 

Jr oxj j^gR V oxj J Jr 


■—(a; y) 


dy. 


(27) 


(28) 


(29) 


Proposition 3.3. For each j = .. .n, 


sup 

rc^K”' 


(a^ y) 


dy < 00 . 


Proof. We first assume (and prove later) that we can differentiate under the integral 
sign, that is 




dQ 


n.-l dx 


-(a;,y) dy2 .. .dy„. 


(30) 


By the Harish-Chandra formula (22), Qi{x^y) can be written as 


Qi(a;,y) = ( 27 r) "/^c„ 
= (2^)-/2c„ 



- ^Tr {Y - UXU^f'^ dU 

- ^Tr {Dy-UDX^f) dC/, 


where D^, Dy are diagonal matrices with the eigenvalues of X and Y as its entries 
respectively. The second equality follows from the first due to the invariance of Haar 
measure on U{n). Observe that by the cyclic property of the trace and the fact that U 
is unitary, Tr {Dy — UD^U^)'^ = Tr {U^Dyll — D^)^. Expanding the trace inside the 
exponential we have 

Tr {Dy - UDX'^f = Tt:DI + Tr DI - 2Tr D^U'^DyU. 

Therefore, 

^{x,y) = c; / X{yf{{U^DyU)yy - Xy) exp ( - ^Tr {D, - DyUf) dU, (31) 

Ju(n) \ I / 
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where = (27r) For a Hermitian matrix H, one can check that Tr = 

Eti hi + 2 Ez<, and therefore Tr [D, - DyUf = YJl=i (^^ " {U^DyU)uf 

+ Then, 


dQi 


dx 


{x,y) 


i=l i<j 


(32) 


where C = 2 sup,j,gR xe = y^2/e. Hence, 


/> (* (* 1 ^ ^ 

/ dyi < Cc), / / A(y)^ exp (--Tr ([/■''dC/]^d?/i. 

Jb ciajj jRn \ 4 ' 

We can make a standard change of variables to the space of n x n Hermitian matrices 
'Hin) by the rule dF = dydU where Z„ = and dF is the product 

of Lebesgue measures Ok^ ^Vij ni<j ‘^Vji- The right hand side of the previous display 
is then equal to 


2-»/2^-n^/2 


5-Tr(Y-DAV4 


JH{n) 


J|g-(yii-:j:i)V4 JJ'g-(y?,-+y|.)/2 < 2’^V2^ 


i<j 


It remains to justify the swapping of the derivative and the integral in (30) and (31). 
For this we shall use the following result from [Bil95, Theorem 16.8]. 


Proposition 3.4. Let (F,/i) he a measure space. Suppose that f{x,y) is a continuous 
and integrahle function of y for each x G I, where I can be taken to be K and that for 
each y GY, ^ {x,y) exists. If for each x* there exists a function g{x*,y) integrahle in 
y such that |§“(a^)2/)| < 9{^*jy) for all y and all x in some neighbourhood of x*, then 

fix,y) M(d2/) = y) y-i^y)- 

Thus, we need to show that Qiix, y) satisfies the hypothesis of the above proposition. 
Since the function x i-A p'^it, x, y)/A{x) A{y) is smooth on R", the same is true for 
Qtix, y) so it remains to find a dominating function g. 

Firstly, for (31), one can apply Proposition 3.4 with g equal to a constant since 
g-Tr (D„-t/DxC/t )=^/2 y ^ and U{n) is compact. For (30), consider the interval \x* — 
h, X* + h] around a fixed point x* gR where h > 0 . Then for Xj G [x* — h,x* + h], we 
have 

— g-y|/ 2 g-a:f/ 2 g 2 ;jyj ^ g-i/J/ 2 g(a;*-|-Ii) I 1/3 | _ ^-(yj-{x’’j+h)f /2 ^(x* + hf I2 ^ 


Therefore, for such Xj, we have by the bounds (32) and (23) that 

dQi, 

■{x,yj \ 

lu{n) 


dx. 


Cc'^ [ A{yf exp ( - ^Tr {U^ DyU - D^f) dU 
Ju(n) \ 4 : J 


<CS^A{yf\{ 


g-(yi-a:i) /4g-(yj-(a;*-|-Ii))^/4g(rc*-|-?i)^/4 


=: 9ix*,y), 
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and g is integrable over with respect to y 2 ,---,yn due to the Gaussian factor. 

Considering j/i, Xi, i ^ j fixed and applying Proposition 3.4 with the above g proves 
(30) and hence completes the proof. □ 


Proposition 3.5. For all j = 1,... ,n 


sup 

(a:,y)GR’*xR 


■—(a; y) 


< oo. 


To prove this we shall use the following formula for the one point correlation function 
K. For 1 < TV < n it was shown in [JohOlb, Proposition 2.3] that the TV-point correlation 
function of Qt is given by a determinant: 

Qt{x,y) dyAT+i ...Ayn = det 

where 


n! 


(n — N)\ J^n-N 


Kt{x,u,v) 


1 

(27ri)^t 



die 




W — Z 


n 

t=i 


W — Xj 
z-Xj 


(33) 


where 7 is a closed contour around the Xt^s and Tl : t ^ L + it, t € M. with L G R large 
enough so that 7 and Pi do not intersect. Then K{x, y) is simply Ki{x, y, y). It is 

sometimes convenient to use the following alternate expression for Kt, see the equation 
below (2.18) in [JohOlb]: 


Kt{x,u,v) = - 




dz din 
l-f Jtl 


n 


W — Z Z — Xj 


n , X 

. , . , Xj(w — z) 

(w-\-z)(w — z)-\-uz — vw-\-1 y - - —r--- 

^_^{w-Xj){z-Xj) 


(34) 


with the same contours as before. Observe that the integral formulas (33) and (34) make 
clear the symmetry of Kt with respect to the ordering oi xi,... ,Xn and that there are 
no issues if any of the Xi's coincide. 


Lemma 3.6. For all x € R" and y G R 


dK, ^ I f dz f 

dxj ^ n 2m Jro 


dw e- 0 - 0 V 2 e(™-y )"/2 
27rz (z — Xj)^ 


n 


W — Xi 
Z- Xi' 


(35) 


Proof. Since 


d n 

-p-r W — Xi W — Z -p-r W — Xi 

dxj z — Xi (z — XiV z — Xi' 

■J 1=1 i/j 

the derivative with respect to Xj of the integrand in the formula for K{x,y) is equal to 


f{x-,z,w) 


1 e-(^-i')V2e(t»-y)V2 
n (z-Xj)^ 


n 


W — Xi 

z- Xi 


The rest of the proof is devoted to justifying the exchange of integral and derivative. 
Consider a bounded set B in the complex plane and let x = (xi,..., x„) with the x^’s 
all lie on the real line in B. Let 7 be a closed contour containing B and therefore also 
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contains a;, then there exist constants d > 0, C > 0 such that for all z G 7, \z — Xi\> d 
for all i and \z\ < C. Moreover, 


W — Xi 


z- Xi 



z — Xi w — z 



w — z 

= 

1 + 


Z- Xi 


Z — Xi Z — Xi 

Therefore, for all a; G B there is a constant such that 


< 1 + 


\w\ + z 


( 36 ) 


\f{x-,z,w)\ < -^snp\e !/)"/ 2 ||g(«« yf/'^\(^(d + C)^ ^ 
=: g{z,w). 

The function g is integrable along the contours 7 and T l. Indeed, 


\w 


n—1\ 




f dz f die sup|e ^^Kd + C)"' ^ 

d-y Jrr. ze7 

(d + C)”-isup|e-("-^)'/2| f din 
267 dr„ 


/Fi Z67 

&„length(7) 


p(t"-y)^/2| 


d"+i 


where in the last line we have shifted the contour Tl to T^ : t ^ y + it hy Cauchy’s 
theorem. The integral with respect to w is just a Gaussian integral and integrates to a 
constant. The other term is treated in a similar fashion but the dm integral is instead 
equal to 

f dw = f I2/+ dt < 00, 

«/ r y «/ ]R 

for each fixed y G M. Thus, by Proposition 3 . 4 , we can differentiate under the integral 
to see that the derivative of K{x, y) is given by 


dK 


1 


— (a:,y) = - / — 
oxi n ' 


dz 


2Tri Stti 


dm e y)^/2e(t" y)^/2 ^ y; _ 2 -. 


{z-XjY 




z - Xi 


Finally, by Cauchy’s theorem we can shift the contour T^ so that L = 0 since there is 
no longer a pole at z = m. □ 


Proof of Proposition 3.5. It is clear from the contour integral (35) that ^-(x, y) is trans¬ 
lation invariant in the sense that §^ix + hi, y + h) = ^^{x,y) for all d G M. Hence, 
sup( 2 ,_j,)gjjnxR is equivalent to sup 3 ,g][jn ^^(x,0) so we only need to bound the 

latter. Fix a constant d > 0. By Cauchy’s theorem, we can take 7 to be the closed 
(rectangular) contour around xi,... ,Xn composed of four parts jt, 7 b, 7 r and 7 ;, where 
jt ■ u ^ —u + di, u G [—i?, B], 76 : u ^ u — di, u G [—B, B], 7 ,. : v ^ R + vi, v G [—d, d], 
and 7 / : u —> —B — vi, v G [—d, d]. B := R{x) is chosen so that the minimum distance 
between the contour 7 and the Xi’s is at least d. We shall consider each parts of the 
contour separately. Denote the integral along the contour 7 ^ by I{yt) and likewise for 
the others. 

Since \z — Xi\ > d for all i and z Gy, we have by (36) that 


On 7r, |z| = |B -f vi\ = {R? + < (B^ -I- d^)^/^ and 

|g-zV2| ^ |g-(R^+2yR«-uO/2| < ^-R^/2^d^/2^ 



W — Xi 


Z - Xi 


< 
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Therefore 


+ ^ dt 

on—2 

length( 7 ,)(d + (i ?2 + d 2 )i/ 2 )"-ie-«V 2 gdV 2 
length(7j.)e" 


(27r)3/2dn+i 

C„2"-2 

(27r)3/2d"+i 

where length( 7 r) = 2d and 


,-RV2„d /2 


a = 


V^Jr \2-/2(i(n-l 


if n odd 
— 1 ))! if n even 


n> 2. 


(37) 


(38) 


Due to the exponential term e ~^we see that the two terms on the right hand side of 
(37) vanishes as i? —> oo and hence the same is true for I{"/r)- By symmetry, the same 
argument shows that /(yz) also vanishes as i? —> oo. Thus, we can deform the contour 
7 to the two horizontal lines, 7 + : u —>■ —u + di and 7 _ : u —> u — di, u £ R.. 

On 7 +, \z\ = (u^ + d^)^/2 amj |g-z /2| _ |g-(-ii+di) / 2 | ^ ^-u / 2 ^d / 2 ^ Hence, in a 
similar fashion as above, we have 


S(7+)l < 


-yn—2 


»dV2. 


1 


27rd"+i 

^ 9n-2 

d /2 


[ (rf+(^2^^2)l/2^"-lg-«V2 du 

Jr 


27rd"+3 
2”-2 


27rd"+i 


•\/^ J B 




where 


c: = 


1 


< 


< 


V^. 

2«-2 

2«-2 


(d+(y2+d2)l/2^-lg-uV2 
(d"-l + (u2+d2)("-l)/2)e-“V2 d^ 


d«-lg-uV2 20-3)/2(y»-l ^n-l^g-uV2 d^ 
2 n—22(71— 3)/2 


= 2”"2d”"i(l + + 






du, 


and the integral on the last line is equal to zero if n is even and equal to (n — 2 )!! if n 
is odd. By symmetry, the same bound applies for /(y-) and hence we have shown that 
there exists a constant C depending only on n and d and is independent of x such that 


sup 


dK, 


< sup (|/(7-i)| + |/(7-)|) < C, 
£ceR" 


as required. 

We now turn our attention to showing (27). Observe that 

/ l^(a;, 2 /)^ dy < sup7ir(a;, y) / K{x,y) dy = nlsup K{x,y), 
Jr y^R Jr y^R 


□ 


( 39 ) 


















CONTINUITY AND STRICT POSITIVITY OF THE MULTI-LAYER EXTENSION OF THE SHE 19 


since Qi{x, y) dy = 1 for all x. So it suffices to show that sup^ j, K(x, y) is bounded 
or equivalently by the translation invariance of K which follows from the translation 
invariance of Qt that sup^jg^n, K{x,0) is bounded. 


Lemma 3.7. 


sup K(x,y) = sup K{x,d) < oo. 

a;GR" xGR" 


Proof. It is convenient to use the contour integral formula (34) instead. Notice that 
there is no longer a pole at in = z and so we can deform the contour so that L = 0. 
Let 7 be the contour in the proof of Proposition 3.5 comprising of four parts, 7 ^, -ji, 
7 t and 7 ^. It can be shown in the same manner as in the proof of Proposition 3.5 that 
the contributions from 7 ^ and 7 ; vanishes at infinity in the direction of the real axis 
and so we can deform the contour 7 to the two horizontal lines, 7 - 1 - : it —> —it + di and 
7 _ : u —>■ It — di, u G R for a fixed d > 0. We then have 

n 

W — Xd 


'''<"•'’> = -72^5/ '>=/ + =)n 7 r 

[ZlTl) JYq • ^ 


Xi 


dz 


die Y[ 


n 

W — Xj 


Xi 


(27ri)2 4 ^u7- Jro ^ “ ^3) 


—: Ii +12- 


(40) 


Denote the contribution from 7 + by j = I, 2 and likewise for 7 _. Note that on 

7 -i_, |z|2 = {vf +d^), — \ > d for all j and z G 7 - 1 -. Hence, 

by (36) we have in a similar manner to the proof of Proposition 3.5 that 

f p 5 / |t| -|- (^2 


|/i(7+)l < ^^dz^dt e-“^/2e-*^/2(|t| + (it^ + d2)V2)(^i + 


■d2)l/2 


(41) 


for some constant Cd,n- By symmetry /i( 7 -) is bounded by the same constant. 
It remains to bound 12- Observe that 


n 


W — Xj 


Xk 


z-Xj ^^{w- Xk){z - Xk) 




W — Xj 


1 


z- Xj ^Z-Xk 
k—1 


^ I . 

< ^ 1 + 


w 


Thus in the same way as above, both |/ 2 ( 7 -|-)| and |/ 2 ( 7 _)| are bounded by some con¬ 
stant „. Combining this with (40) and (41) shows that there exists a constant C 
independent of x and depending only on n and d such that 


sup K{x, 0) < C, 


which completes the proof. 


□ 


Proof of Theorem 3.1(a). Lemma 3.7, Proposition 3.3 Proposition 3.5 and (29), (39) 
together imply that (27) and (28) are bounded. This in turn shows that the assumption 
of Lemma 3.2 is satisfied and the result follows. □ 


Lemma 3.8. For all t > 0 and x G 
on n such that 


there exists a constant C 4 > 0 depending only 

f Kt{x,y)‘^ dy < 

Jr 
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Proof. By the scaling property of Qt and a change of variables 

[ Kt{x,yf dy = r^^'^ [ Ki{xt~~^/'^,y'f dy'. 

By Lemma 3.7 and (39), the latter integral for each fixed n is bounded uniformly in x 
which gives the desired result. □ 

Proof of Theorem 3.1(b). Let t = u + h where h > 0, then we need to estimate 

/ / {Ks+hix,y) - Ks{x,y)Y dyds. (42) 

Jo Jm. 

Assume for now that one can differentiate under the integral in formula (33). The time 
derivative of Kt is then equal to 


d 

—Kr{x,y) 


dz f dm 1 ^ 


w — Xi e 


{w-yp /It -{z-yf I2r 


'Tl 


2Tri 


i=i 


2 — Xi 


r'^(w — z) 


1 + 


(m - j/)2 (z - yf 


dz f dw TT 

Otto /^, 9 ' 7 r 7 J-A 


w — x) 6“ 


r-v/r A/ 2Tri /p/ 27ri z 
' r j=l 


{w — z) 


2r 


_ w z 

^ + T“T 


where x) = (xj — y)!\pr and 7 ', T'j^ are the contours 7 , Li translated by y and scaled 
by Ijyjr. We can rewrite the derivative as 

^Kr{x,y) = - ^-j=Ki{x',d) - \=I{x') 

or Tyjr Vyjr 


where 


I{x') := 


" dz 
v/ 27r* 


— TT 

27ri z — xl 
f=i J 


- X[ 


{w — z) 


(w^ - z^) 


w — X 


'' dz r dw ^ 

Y 2711 7p/^ 2717 z — X 


/\w + z), 


where in the last line we have shifted the contour L^ so that L = 0 as there is no longer 
a pole at m = z. 

Note that \{w + z)e“ /'^\ is uniformly bounded on the chosen contours and as 

in the proof of Proposition 3.5, we can deform 7 to the two horizontal contours 7 -|_ and 
7 _. Thus, there exists a C := C"(n, d) such that 




dtt; 


— Xa 


|e“ 


Essentially the same calculation as for Proposition 3.5 shows that sup 3 ,gRn |/(a;)| < 00 
and together with an application of Lemma 3.7 gives 


Krix,y) 


<J:^{\K{x',0)\ + \I{x')\)< ^ 


v- ’-U ' H V- /I/ 

for some constant C independent of x. 

Now, rewrite the integrand in (42) as \Ks+h{x,y) — Ks{x,y)\‘^ = \Ks+h{x,y) — 
Ks{x,y)\‘^~°'\Ks+h{x,y) — Ks{x,y)\°' for a < 1/2. We estimate the latter factor by 

< C'“(syi)-“|/i|^ 


fs+n Q 

J -^^r{x,y)dr 
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For the other term we have by time scaling 
\Ks+h{x,y) - Ks{x,y)\'^~°‘ 

< \Ks+h{x,y) + K,{x,y)\‘^~°^ 

I _ _ 1 

= . K^x/Vs + h,y/Vs + h)-\ — ^K{xly/s,yjy/s) 

y/s + h vs 

X {Ks+hix, y) + Ks{x, y)) 

^ 2i-“s(“-i)/2( sup K{x, 0 )Y °‘{Ks+h{x,y) + Ks{x,y)). 
a;GR" 

Therefore, foi u <t <T, the right hand side of (42) is bounded above by (the constant 
C := C{n, a) may change from line to line) 

C\h\^ f s-3“/2s(“-i)/2 f Ks+h{x,y)+K,{x,y)dyds 
Jo Jr 

nU 

= C|/l|“ / s-(“+l/2) ds 

Jo 

< cr5-“|/i|“, 

since a + < 1. 

It remains to justify the differentiation under the integral sign in -^Kr{x,y). We 
once again appeal to Proposition 3.4, which means finding a dominating function g for 
the derivative. Let /(r; w, z) denote the integrand in (33) ioi u = v = y for fixed y S R 
and X € R" (note that we have suppresed the dependency on x and y in the notation), 
then differentiating with respect to r we have 

d , 1 -^W-Xj 

—f{r;w,z) = — M - - - 

or z — Xj w — z 

j=i J 

n 

+ rr + z-2y) 

2j-i J.X z — Xi 
j=l J 

=: Ii +12 

Let zf- = sup^g..^, Re(z) and zl = sup^^g^ Im(z). Fix r* > 0 then for all r G [r*/2,2r*] 
and z G 7 , we have that 

|g-V-y)V2r| < 


and for all w G F^ : 1 1 -^ L + it, we have 

|g(>"-!/)V2»'| < l^r, jr, ^ 

Hence, for all r G [r*/2, 2r»] 

w — X 


Ii4l<4rr^n 


rf \w - z 




Z — Xi 


g-V4r. g-tV4r. g(L+yf/r 


. n 

* i=i 

Let g(rt,;w,z) be the sum of the upper bounds of |Ii| and II 2 I, then it can be shown in 
a similar fashion as in Lemma 3.6 that ^(r, ;w,z) is integrable on the contours 7 and 
F L and so an application of Proposition 3.4 completes the argument. 


2?/l 


/4:r, (zlf /r„ /ir, JL+yf /r 
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Finally, by Lemma 3.8 we have 



Kt-s{x,yY Ayds < Ca f {t — s) As < 2CA\t — . 

J U 


This completes the whole proof of the theorem. 


□ 


4. Existence, Uniqueness and Moment Estimates 

4.1. Bounded Initial Data. We now prove the existence, uniqueness and moment 
estimates part of Theorem 1.3(a). The proof of continuity will be delayed to Section 5. 
In the sequel constants will generally be denoted by c, C or K and possibly adorned 
with primes, tildes, subscripts or superscripts. They may differ from line to line and 
their dependence if any will always be specified. However, Q, 1 < i < 4 will always 
mean the constants in Theorem 3.1 and Lemma 3.8. T > 0 will always denote the finite 
time horizon. 

Proof of existence, uniqueness and moment estimates of Theorem 1.3(a). The proof is 
by a Picard iteration argument. Throughout the proof, we fix an arbitrary integer 
p>2. For {t,y) G (0, oo) x K" define m^{t,y) := Jn{t,y) where was defined in (15) 
and for fc > 1, let 

m’'{t,y) = m°{t,y) + An [ [ Qt-s{,y,y')m'"~^{s,y') Ay'., W{As,Ay[) 

Jo Jr'^ 

=: m°{t,y) + l'^(t,y). (43) 

We first show that each of the stochastic integrals above are well defined, that is for all 
(t, y) £ (0, oo) X R", the random field (/fc(s, x), (s, x) £ (0, t) x R) defined by /fe(s, y{) := 
J^n-i Qt-s{y,y')m’^{s,y') Ay( is in ^2 for all fc > 0. 

Fix {t,y) £ (0,oo) X R"- and consider fo(.s,yf) = /g„_i Qt-s{y,y')'mP{s,y') Ay(. We 
need to show that mP satisfies the three assumptions of Proposition 2.5. Since the initial 
data g is .^o-measurable, mP is adapted to the filtration By assumption on g, 

suPygRn ||5'(j/)|lp < Kp,g < 00 and hence by Minkowski’s integral inequality 

\\m°{s,y)\\p < —^ f \\9{y')\\pQt{y,y') Ay' 
ni Jg^n 

< sup ||5(2/)llp^ [ Qt{y,y')Ay' 

< Kp,g. (44) 

Therefore, ||to°(s, y)\\p is bounded by K^ g uniformly for (s, y) £ [0, 00 ) x R". By Lemma 
5.2 below, {s,y') 1 —>■ mP{s,y') is continuous in Lp(yi.) on (0,t) x R" and so Proposition 
2.5 implies that /o £ 3^2 and 

[ f Qt-siy,y')m°{s,y') Ay',W{As,Ay'0, 

Jo Jr^ 

is a well-defined Walsh integral. Consequently, the random field {w}{t,y) = m°{t,y) + 
I^{t,y),{t,y) £ (0,oo) X R”) is well defined. 

We wish to show that the sequence {m'^{t,y)}k>o is Cauchy in L^(f2). To this end, 
let dk{t,y) := ||m*’+^(t,y) — m^{t,y)\\p. By Lemma 2.2, Lemma 3.8 and (44), we have 
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for all {t,y) G (0,oo) x R", 

do{t,yf<Alcl [ [ ( [ Qt-siy,y')\\m°{s,y')\\pdyi'] dj/ids 

do JR VdR"-1 / 

< 2Kl^C^AlclVt 


where (74 is the constant in Lemma 3.8 and r(3/2) = 12 . 

Now assume that for all 0 < Z < fc, y), {t, y) G (0, oo) x R”) is well defined and 

satisfies 

(i) m} is adapted, 

(ii) (s, y) !->■ w}{s,y) is L^(r2)-continuous on (0,t) x R" for all t > 0, 

(iii) for all (t, y) G (0, oo) x R" and 0 < Z < fc — 1 


di{t,yf < Klg{CiAlcly/^)^+^ 


ia+i)/2 

r(^ + i)' 


We want to show that the same is true for and dk- Let (Z, y) G (0,oo) x R". 

Observe that m^{t, y) = wP{t,y) + Yl\=i v) — y), and so to bound the pth 

moments of it sufhces to bound each of the di’s, 0 < I < k — 1. Indeed, by property 
(iii) and (44), we have 


\m^{t,y)\\l < 2||m°(Z,y)||2 +^2'cZ;_i(Z,y)^ 




^1/2 




(45) 


which shows that sup(g_j,)g[g ||w^(s,y )||2 < oo. This and the induction hypothesis 
shows that satisfies all three assumptions of Proposition 2.5 and so fk G ^2 and 

l''^\t,y) = f f Qt-s{y,y')m^(s,y') dy' W(ds,dy'i), 
do dR"-l 

is a well-defined Walsh integral for all (Z, y) G (0, 00 ) x R". Moreover, it is adapted and 
so 771^+^ = m° -I- is also adapted. We need to check the L^(f2)-continuity of 1^+^. 
Fix a < 1/2, then for all 0 < r < u < Z and y, 0 G R” 

||/'=+i(u,y)-/'=+i(r,z)||i 

< / f ( f {Qu-siy^y’) - Qr-.s{z,y'))\\m'^{s,y'\\ 2 dy’^'] dy'^ds 

do dR \dR"-l / 

+ 2An / [ ( [ Qu-s{,y,y’)\\m^is,y')h dyi') dy/ds 

Jr dR \dR"-i / 

< 2 Al(Ci+C 2 + C^) sup \\'m!‘{s,y')\\l{\y - z\ + \u-r\°‘), 

(s,y')G[0,u]xR" 


by Theorem 3.1 which proves the L^(r2)-continuity of on (0,Z) x R". 





CONTINUITY AND STRICT POSITIVITY OF THE MULTI-LAYER EXTENSION OF THE SHE 24 


For the bound on dk, we use Lemmata 2.2 and 3.8 and the induction hypothesis to 
obtain 


dk{t,yf < 

= Kl(C,AlclV^f+^- 


, fe /2 


/o r(| + i) 

fik+l)/2 




r(^ + i)’ 

where we have used the Euler Beta integral [OLBCIO, equation 5.12.1]: 


du = 


r(a)r(6) 
r(a + b) 


, a,b > 0, 


(46) 


(47) 


and the fact that r(l/2) = y/ir to evaluate the time integral. It follows that the bound 
(45) holds with k replaced with fc + 1 and that sup^^ ||to^+^(s, y )||2 < oo. 

Hence, satisfies all the assumptions of Proposition 2.5 and therefore fk+i G ^ 2 - 

Thus, by induction we conclude that for all integers k, the random field {rn^{t,y) = 
wP{t,y) + l'^{t,y),{t,y) S (0, oo) x R") is well defined and satisfies properties (i), (ii) 
and (iii) listed above. 

We now show that the sequence y)}k>Q is Cauchy in L^iyi). This follows from 

the fact that for any T > 0 


OO 

sup ^dfe(t, y)<oo, 

(t.ylGlO.TlxR" 


which is a consequence of property (iii), the ratio test and the following asymptotic: 

~ as z ^ oo, see [OLBCIO, equation 5.11.12]. We conclude that there 

exist a random field which we denote by Mn{t,y) such that m^{t,y) — >■ M„(t,y) as 
A: —>■ 00 in L'P{yi) and almost surely for a subsequence uniformly in y G R" and t G [O.T]. 

Since each is adapted, M„ is also adapted. The L^(H)-continuity of M„ is inher¬ 
ited from that of since the convergence is uniform on [0, T] x R” for all T > 0. Now 
take fc —>■ oo on both sides of (45). By [CD14b, Proposition 2.2], we know that for all 
a; > 0 


(l + erf(a:)) = 


r.k-1 


^r(^) 


(48) 


Using this with x = 2C'47l^Cp-yTt^/^ gives the bound (16) in the statement of the 
theorem. Thus, by Proposition 2.5, for all {t,y) G (0, oo) x R” the random field / 
defined by /(s,yi) = Qt-siy,y')M„{s,y') dyi for {s,y[) G (0,t) x R is in and 
the stochastic integral 


In{t,y) 



Qt-s{y,y'M„{s,y') dy(. lU(ds,dyi), 


is well defined. 

It remains to show that the limit Mn{t,y) solves (15). Fix {t,y) G (0,oo) x R". By 
definition, m^{t, y) = wP{t, y)+I^{t, y) where the left hand side converges in L^(r2) and 
almost surely for a subsequence to Mn(t,y). For the right hand side we have by the 
uniform convergence U'{Vl) of that 


Wl'^ii^y) - In{t,y)\\l<‘2.ViAl^cl sup ||TO'=(s,y') - ^n(S:J/')llp 

(s,?/^)G[0,i]x]R”' 

—>■0 as fc —>• 00 . 
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Therefore, we have LP{VI) convergence of to In{t,y) and hence almost sure 

convergence for a subsequence to the same limit. The limit of both sides of m^{t,y) = 
mP{t, y) y) must be equal almost surely and so we have shown that for all {t, y) G 

(0, oo) X R", Mn{t,y) satisfies (15) almost surely. This proves existence. 

For uniqueness, suppose that M{t,y) and N{t,y) are both solutions to (15) with the 
same initial data g and let d{t,y) = \\M{t,y) — N{t,y)\\p then by a similar calculation 
as for existence we have 


d{t,y)'^ < 


sup 

(s,y)e[ 0 ,t]xM^ 


d(s, yf {C4Alcl^/^) 


^n/2 


(49) 


which converges to 0 as n —> oo since the expression on the right hand side is summable 
in n. Therefore, d = 0 and so for all {t,y), M{t,y) = N{t,y) almost surely i.e. M and 
N are versions of each other. This proves uniqueness. □ 


4.2. Delta Initial Data. 


Proof of existence, uniqueness and moment estimates of Theorem 1.3(b). Fix an integer 
p > 2. We first show that if solutions to (13) exists then it must be unique. Suppose 
M{t,x,y) and N{t,x,y) are two solutions to (13) and let d{t,x,y) = \\M{t,x,y) — 
N{t,x,y)\\p. By linearity of the equation (13), M{t,x,y) — N{t,x,y) is a solution to 
(15) with zero initial condition i.e. M{t,x,y) — N{t,x,y) = M-l{t,y) with g = 0. Then 
by (16), sup^, ygRn (i(t, x, y)^ is a bounded function of t S [O,?^] for any T > 0. The 
bound (49) applies to d{t,x,y)‘^ which shows that M{t,x,y) = N{t,x,y) almost surely 
for all {t,x,y). This proves uniqueness. 

We now prove existence. We shall show that Mn{t,x,y) defined by (11) satisfies 
equation (13) for all (t,x,y) G (0,oo) x K" x K". Recall that Mn(t,x,y) is well defined 
on the boundary of the Weyl chamber and it is symmetric under permutations of both 
its space variables, hence we can extend it to a function on R." x R". Similarly we 
also extend Qt-s{x,y) to the whole of R" x R”. Substituting the chaos expansion of 
Mn into the stochastic integral term of (13), using the expression for the correlation 
function Rk (20) and the stochastic Fubini’s theorem [Kho09, Theorem 5.30], we have 
bearing in mind that we can interchange the summation and integral because the series 
is convergent in L^(r2) that 


An f f Qt-si{y,y^)Mn{si,x,y^) dylW{dsi,dyl) 
Jo 


— An 


Qt-sAy,y^)p*ni.si,x,y^) , i 


A(x)A(yi) 


dyi lT(dsi,dyi) 


f P*nit-si,y,y^) ^ 


fc+l 


/o 


A(a;)A(y) 




Wp*n{si-i-Si,f \y*) 


fc + 1 


X Pn{sk+I,y''""^,x) P[ dy* W®'‘(ds,dy) dyl VF(dsi,dyi^) 


2=2 


P*n{t,x,y) 

A(a;)A(y) Jq 


[ [ Ri{si,yl-,t,x,y)W{dsi,dyl) 
Jo Jr 


oo p p fc+l 

+ X! / / A^+^pl(t-si,y,y^)Y[pl{s,-i- Si,y^-\y^) 


k+1 


X pUsk+1 , y"+', x) n dy: (ds, dy) 


2=1 
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P*u{t,x,y) f r 


Rk{s,y,t,x,y) dy), 


where the last equality follows by a relabelling of the indices. Thus, the right hand side 
of (13) after the substitution is equal to 


P*u{t,x,y) 

A{x)A{y) 



Rk{s,y;t,x,y) lT®'=(ds,dy) 


which is the definition of Mn{t, x,y) as required. 

It remains to estimate the pth moments of M„(t, x, y). The approach is to construct an 
approximating sequence to M„ and estimate the moments of each term of the sequence 
and take limits. The natural candidate for the approximating sequence is the following: 
for each {t,x,y) G (0,oo) x R” x R”, let m°{t,x,y) := Jn{t,x,y) where Jn was defined 
in (13) and for fc > 1 define 


{t,x,y) =m°{t,x,y)(l + '^ [ [ Ri{s,y-t, x,y) W®\ds,dy) 


In other words, m^{t, x, y) is the fcth partial sum of the chaos expansion for Mn{t, x, y). 
Let dk-iit, X, y) := x, y) — x, y) ioi k > 1, then clearly 


By Lemma 2.3 


dfc_i(t,x,y) = m°(t,x,y) / [ i?fe(s, y; t, a:, j/) lT®''(ds, dy). 

JAk{t) Jr'‘ 

2.3 

\\dk-i{t,x,y)\\l<cl'‘m°{t,x,yf [ [ Rk{s,y-t,x,yf dyds. 

aAfe(t) 


It is easy to see that 

k 

and so by (50), we have 

k 

X, y)\\l < x, yf + ^ 2'||di_i(t, x, y)\\l 


(50) 


< 


2m°{t,x,yf(l + ^{2ciy[ [ i?i(s,y;t,x,j/)^ dydsV 

V JAiit) Jr‘ J 


Each term in the sum above is equal to (2Cp)^Ef^y^j [(— /l\ by Lemma 

2.6 where X = (X^,..., X"), Y = (E^,..., F") are independent copies of a collection 
of n non-intersecting Brownian bridges which start at x in time 0 and end at y in time 
t. Letting /c —>■ oo we have for all (t, x, y) £ (0, oo) x R" x R" 


lim ||m''(t,a;,j/)||^ < 2TO°(t,x,y)^E 

k—^oo ^ 


X,Y 

x,y;t 


exp (2cl J2 Lt{X^ - Y^)) 
* j'=i 


(51) 


For each t > 0, Lemma 2.7 shows that the right hand side of the above is bounded 
uniformly in a;, y G R" for any p > 2. By Cauchy-Schwarz inequality 

X, y)-m!^' {t, x, y)\\P < x, y)-m^' {t, x, y)|| 2 ||m'=(t, x, y)-m^'(t, x, y)|l 2 (^’_i), 

which converges to 0 as k, k' —>■ oo by the L^(n) convergence of and the moment 
bound (51). Therefore, m^{t,x,y) also converges to Mn{t,x,y) in LP(n) and we can 
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replace the left hand side of (51) with \\Mn{t, x,y)\\p. This completes the proof of 
existence, uniqueness and moment estimates. □ 

5. Continuity 

We shall use the following version of Kolmogorov’s continuity criterion which is due 
to Chen and Dalang, see [CD14a, Proposition 4.2]. 

Theorem 5.1. Consider a random field {f{t,y) : {t,y) G M+ x R’’*}. Suppose there are 
constants oq, ... ,ad € (0,1] such that for all p > 2{d + 1) and all M > \, there exist a 
constant C := C{p, M) depending on p and M such that 

\\f{t,y) - f{s,x)\\p <c(\t- sl““ + X! “ 2;^]“’ 

i=l 

for all {t,y) and {s,x) in [1/M,M] x Then f has a modification which is 

locally Holder continuous on (0, oo) x R'’* with indices {fiag, ..., fictd) for all fi G (0,1). 

5.1. Bounded Initial Data. We now prove the Holder continuity of the solution to 
(15) by verifying the assumptions of Kolmogorov’s continuity criterion. We first esti¬ 
mate the increments of Jn{t,y) = ^ g{y')Qt{y,y') dy' where g satisfies the bound 

SUPygRn \\g{y)\\p < Kp^g. 

Lemma 5.2. Let M > 1 and p > 2. There exist constants Ki := Ki(M,n,p) > 0, 
7 = 1,2 such that for all t, t' G [i/M, M] and y, y' G R" 

\\Jn{t,y) - Ju{t',y)\\p < Ki\t-t'\, 

and 

\\Jnit,y) - Jn{t,y')\\p < K2\y - y'\. 


Proof. By the assumptions on g and Minkowsky integral inequality, we have 

\\Jn{t,y)-Jn{t',y)\\p< ^ sup \\g{z)\\p [ \Qtiy,z)-Qtfiy',z)\dz. 
n'. zgR" Jr” 

For t > 1/M, Qt has bounded derivatives in both time and space and the result follows 
by a direct calculation. □ 

We now turn our attention to the stochastic integral term /„(<,j/). 

Proposition 5.3. Let M > 1, a < 1/2 andp>2. For all (t, y) and (u, z) G [0, M] x R" 
there exists a constant K := K{a, g, M,n,p) such that 

\\Init,y) - Iniu,z)\\p <K{\t- Jy- 


Proof. We consider the spatial and temporal increment separately. By (16), there is a 
constant C := C{g, M, n,p) such that 


sup l|M„(s,y')llp < C'- 

(s,p')g[ 0,M] xR" 
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Then by Lemma 2.2 and Theorem 3.1(a) 


<AWp / / {Qt-s{y,y') - Qt-s(z,y'))\\Mnis,y')\\p dy'A dy[ds 

Jo Jr / 

<CAlcl f [ f [ Qt-siy,y') - Qt-siz,y') dyi') di/ids 

Jo Jr vJr"-! / 

<C,CAlcl\y-z\. 

For the temporal increment we have two terms (assuming without loss of generality 
that 0 < u < t < M) 

\\In{t,y)-In{u,y)\\l<2l + 2ll, 

where by Theorem 3.1(b), for any a < 1/2 there exists a C 2 such that 


I := 


An 


10 JR" 


{Qt-s{y,y') - Qu-s{y,y'))Mr,{s,y') dy' W{ds,dy[) 


<CAlcl [ [ ( [ Qt-siy,y')-Qu-siy,y')dyi) dy/ds 

Jo Jr \Jr"-i / 


< C 2 CAlcl\t-u\^ 

and a constant C 3 such that 


11 := 


f f Qt-siy,y')Mn,{s,y’) dy’^W{ds,dy[) 
Ju Jmp- 

<C:iCAlcl\t-u\^/^ 

<C'3CJl2c2(2M)5-“|t-u|“. 


□ 


By the subadditivity of the function a; i-A |a;|^, for /3 € (0,1] we have 

/ n \PA n 

\y-yf=\^\yi-yi?] <'^\yi-yf- 

\i = l / Z = 1 

Lemma 5.2 and Proposition 5.3 together shows that for all M > 1, a < 1/2 and 
p > 2, there is a constant C := C{a, g, M,n,p) such that for all (t,y) and {t',y') in 
[1/M,M] X 


||M„(t, y) - M„(t', y')\\p < C { \t - t'^/^ + V |y. - 


Taking p large enough and applying Theorem 5.1 shows that M„ has a version that is 
locally Holder continuous on (0, 00 ) x R" with indices up to 1/4 in time and up to 1/2 
in space. 


5.2. Delta Initial Data. We now turn our attention to Mn{t,x,y). Observe that in 
this case we cannot apply the method used in Proposition 5.3 directly since the yth 
moments of M„(t,x, y) are not bounded uniformly in time, for instance if x = y then 


||M„(t,x,x )||2 > 


Pn{t,X,x) 

A(x)2 


(27rt)-’"/2 


n 

(1+ ^ sgn(a)ne-("^-"''«^'/"‘), 

cr^Sn 1 

cr^id 


A(x)2 
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which converges to infinity as t J, 0. However, for any t > 0 fixed we have by (51) and 
Lemma 2.7 that there is a constant C := C{n,p) such that 

2 


\\Mn{t,X, 


\l<2 


P*n(t,x,y) 

A{x)A{y) 


E; 


X,Y 

'x,y,t 


exp (2c^ ^ Lt{X^ - 


< ct- 


uniformly for x, y £ R". Thus, for all positive times, M„ belongs to the class of initial 
data in Theorem 1.3(a). It is clear that at any given time we can restart the equation 
taking the current solution as the new initial data. More precisely, let t > 0 and consider 
the shifted white noise W'^{s,y) := W{t + s,y). Define M^{t,x,y) := Mn{T + t,x,y) 
then it is easy to check by using the semigroup property of Qt that satisfies the 
integral equation 

Mnit,x,y) = ^ [ Mr,{T,x,y')Qt{y,y') dy' 
n\ Jjjn 

+ An f f Qt-s{y,y’)M^{s,x,y') dy'^W'^{ds,dy[). 

Jo Jr" 

In other words, is the solution to (15) driven by the shifted noise with initial 
condition M^(0,a;,j/) = Mn{T,x,y). Now define 


Mnit,x,y) 


Mn{t,x,y) if 0 < t < T, 

M^{t-T,x,y) iit>T. 


Clearly, A4n{t,x,y) solves (13) and by uniqueness, is a modification of the chaos 
series (11). Let M > 1, a < 1/2 and p > 2 then since sup^^ ^ggn \\Mn{T,x,y)\\p < 
oo, Lemma 5.2 and Proposition 5.3 applies to show that there is a constant C := 
C(a, M, n,p, t) such that for all t, t' G [r, M] and y, y' G [—M, M]” and x G R" 

Wl{t,x,y) - ,x,y')\\p < C{\t - (52) 


5.2.1. Continuity in the Initial Condition. We study the continuity of x i—>■ Mn{t, x,y)', 
in fact we show that {t,x,y) i-A M„(t,x,y) is jointly continuous. Recall the chaos 
expansion of Mn(t, x, y): 

Mn{t,x,y) = Jn{t,x,y)(l + '^ [ [ Rfc(s,y';t,x,j/) IT®'=(ds,dy')Y (53) 

V Jr'‘ J 


where for 0 < si < ... < Sfe < t, y = (y(, j/J,..., j/(=) 


Rkis,yC,x,y) 


= At 


Pn(5ua^.l/ )ni=2Pn('Si " s^-l,y^ \ y^)p* (t - Sfc, , y) 

Pl{t,x,y) 


nn<isj 

i=lj=2 


It is easy to see that Jn(t, x, y) = Jn{t, y, x) and from the expression of Rk one can see 
that for all /c > I 

Rk{s,y-t,x,y) = Rk{t - s,y;t,y, x), (54) 

where t-s := {t-Sk, ■ ■ .,t-si), 0 < t-Sk < ... < t-si < t and y := (yf,?/i“\ ... ,y(). 
Therefore, it is reasonable to think that each term in the sum above is symmetric in x 
and y provided one can reverse time in the multiple stochastic integral. This motivates 
the following proposition. 


Proposition 5.4. For all n> 1 and for each fixed y the random fields {Mn{t^x^y), (t^x) 
G (0,oo) X R") and (M„(t, ?/,x); (t, x) G (0, oo) x R”) are equal in distribution. 
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Proof. Fix k>l and {t,x,y) G (0,oo) x R" x R". Recall the time reversed white noise 
W defined by 1R([0, s] x A) = W([t — s,t] x A) for s < t and A G ^b(R). Extend 
Rk{s,z;t,x,y) to a function on x R^) by setting it to be zero for s ^ Afe(t). 

Let Rk be the symmetrisation of Rk given by 

Rk{s,y,t,x,y) = -^ X! x,y), 

■ TreSfc 


where tts = (s,r(i), ..., and likewise for Try. Clearly, we have Rk{s,y\t,x,y) = 

Rk{s,y;t,x,y). Therefore by Lemma 2.1 and (54), (recall the definition of the multiple 
stochastic integral in Section 2.1) 



Rk{s,z;t,x,y) W®'^{ds,dz) = f f Rk{s,z;t,x,y) W®'^(ds,dz) 

JR'' 


[ [ Rk{t - s,z;t,x,y) (ds, dz) 

J[o,t]*' Jr'' 

Rk{s, z; t, y, x) lE®''(ds, dz) 


'[0,*]'' Jr'' 


/Afe(t) Jr*' 


Rk{s,z;t,y,x) lE®''(ds,dz). 


Thus, applying the above to each term of the sum in (53) we see that 


Mn{t,x,y) 


Jn{t,y,x)(l + '^ [ [ Rk{s,y;t,y,x)W®''{ds,dy) 

\ JAfc(t) Jr*' 

Mn{t,y,x), 


for all (t, X, y) G (0, oo) x R" x R” and the result follows. 


□ 


Finally, we return to proving the joint continuity of the solution to (13). We bound 
\\Mnit, X, y)—M.n{t', x', yOlIp by considering the increments in each variables separately. 
Since Mn(t, x, y) = Mf^{t — r, x, y) for t > 2r, we have by Proposition 5.4 and (52) that 
for all M > 1, p > 2 and a < 1/2 there is a constant C := C{a, M, n,p, r) such that for 
all {t,x,y) and {t',x',y') G [2t, M] x [—M, M]” x [—M, M]" 

\\Mn{t,x,y) - Mn{t',x',y')\\p 

<\\Mf;{t-T,x,y) - Ml {t' -T,x,y')\\p + \\Ml (t' - r, y', x) - (t' - r, y', x')Up 

< c{\t - +1^ - + \y- yf/")- 

Since r > 0 is arbitrary, we can take 2r = 1/M and thus we have shown that there 
exists a constant C = C{a, M, n,p) such that for all (t, x, y) and (<', x', y') G [1/M, M] x 
[—M, M]^" the above inequality holds with C in place of C. Finally, using the subad¬ 
ditivity of X !->■ |x|^ for /3 G (0,1] and applying Theorem 5.1 proves the existence of a 
Holder continuous version. This concludes the entire proof of Theorem 1.3. 


6. Strict Positivity 

6.1. A Weak Comparision Principle. Recall that iF„(t,x,y) can be expressed as 
Kn(t,x,y) = det[u{t,Xi,yj)]2j^i where u{t,x,y) is the solution to (2) with initial data 
Sx- Bertini-Cancrini [BC95] proved that u(t,x,y) is the limit in LP(fl) for all p > 2 
of u^(t,x,y) as e —>■ oo, where u^{t,x,y) is the solution to the stochastic heat equation 
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subject to a mollified white noise in place of the space-time white noise. Its solution 
is given by the following Feymann-Kac formula which is well defined for the noise W^: 

u%t,x,y) = pt{x - 

where the expectation is with respect to a Brownian bridge b starting from x at time 
0 and ending in y at time t. By the above Feymann-Kac formula it is then clear that 
for all {t,x,y) G (0,oo) x K x R, with probability 1, u{t,x,y) > 0. Using this and the 
determinant formula for Kn, the authors in [OWll, Proposition 5.5] proved by a path 
switching argument that Kn{t, x,y) >0 almost surely, for all (t, x, y) G (0, oo) x Wn x Wn- 
In fact, a stronger result is true since the above implies that Kn{t,x,y) > 0 for all 
rational points (t, x, y) almost surely. It is well known that it, x, y) i—>■ u{t^ x, y) has 
a jointly continuous version and hence the same is true for Kn as it is just a sum of 
products of the u’s. Therefore, by continuity 

P[Kr„(t, x,y)>0 for alH > 0 and x, ?/ G IF„] = I. 

Since the Vandermonde determinant is non-negative on Wn, we see that the same is 
true for Mn in the interior IU°. By the continuity of Mn proved in the previous section, 
this non-negativity extends to the boundary of the Weyl chamber and by symmetry to 
the whole of R". That is, 

F[Mn{t, x,y) >0 for alH > 0 and x,y £ R”] = I. (55) 

By the linearity of the equation (13), the non-negativity property above is equivalent 
to a weak comparison principle. The next lemma extends this to solutions M^{t,y) of 
equation (15) with initial data g. 

Lemma 6.1 (Weak comparison principle). Let M^(t,y) and Mn{t,y) be the solution to 
(15) with symmetric initial data gi and g 2 respectively. If gi > 52 , then 

F[Mn{t,y) > Mn{t,y) for all t > 0 and y G R"] = 1. 


i?xp( / IU'^(s, 6s)ds 


Proof. By linearity of the equation (15), it suffices to prove the lemma in the case g > 0. 
For {t,y) G [0,oo) x R", define 


Vg{t,y):=-^ [ g{x)Mn{t,x,y)A{xf dx. 
nl Jb" 

A direct calculation shows that Vg satisfies (15) and so by uniqueness Vg{t, y) = Mf^{t, y) 
almost surely for all {t,y) G [0,oo) x R". Now by (55) and the non-negativity of g and 
the Vandermonde determinant it is clear that for all {t,y) G [0,oo) x R", Vg{t,y) > 0 
almost surely. This and the continuity of {t,y) >->• M^{t,y) shows that P[M®(t,y) > 
0 for alH > 0 and y G R"] = 1 as required. □ 


6.2. A Strong Comparison Principle. We now prove a strong comparision principle 
of which Theorem 1.4 is an easy corollary. 
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Theorem 6.2 (Strong comparision principle). 

(a) Let M^{t,y) and M‘^{t,y) he two solutions to (15) with initial data gi and g 2 
respectively where gi and g 2 are as in Theorem 1.3(a). If furthermore gi > g 2 
and gi{y) > g 2 {y) for some y £ R” almost surely, then 

F[Al){t,y) > M^{t,y) for all t > 0 and y G R”] = 1. 

(b) Let Mn(t,x,y) be the solution to (13), then 

F[Mn{t, x,y) > 0 for all t > 0 and x,y £ R"] = 1. 


We begin with a lemma which provides a lower bound for the deterministic term 
Jn{t,y) in (15). 

Lemma 6.3. Let j3 := l3{n) = Pgue[L : <(>i(T) > 0,Vi]/2 > 0 where (fifY) is the ith 
eigenvalue of an n x n matrix Y from the Gaussian Unitary Ensemble (GUE). For all 
h > 0, t > 0, M > 0, there exists an mo '.= mo{h, M, n, t) such that for all m > mo, all 
s £ and x £ Wn, 



Q y')^( — h,h)^iy') dy ^ /51(—li—(^)- 


Proof. Since Dyson Brownian motion is realised as the eigenvalues of Brownian motion 
on the space of n x n Hermitian matrices Tl{n), we have that 

[ (5s(a:,j/)i(-/i,ii)-(y) dy = / Ps(y)i(_^,,,)n((/)(r+ 113,)) dy, 

JW„ JH{n) 

where Ps{A — B) = for A,B £ TL{n) is the transition den¬ 

sity of Brownian motion on the space of Hermitian matrices and (j) : 'H{n) —>■ Wn is such 
that (fix) = y = (yi,..., y„) = (^i(y),..., (fnX)) is the vector of ordered eigenvalues 
of y. Dx is a diagonal matrix with entries x = {xi, ..., x„). Weyl’s eigenvalue inequality 
[Bha97, Theorem IIL2.1] implies that for two Hermitian matrices A, B with eigenvalues 
4>i{A) and (j)i{B), 1 <i <n respectively, the following hold 

(j)i{A + B) < <pi{A) + (j>{B) and (pn{A) + (j)n{B) < <pn{A + B). 

Therefore 


M-h,h)’^{XY + Dx)) = 1{MY + Dcc) > -h}l{MY + Dx) < h} 
> l{(j)n{Y) YXn> -h}l{4>i{Y) +xi < h}, 

and hence 

Qsix,y)l(-h,h)’-iy) dy 


'Wr, 


> 


[ PsiY) l{(/)„(y) > -h- Xn]l{(l>i{Y) <h-xi] AY 

I P,iY) i{My) > ^}i{<*.(r) < di 


/«(« 


n 

piiY) ( 


— h — Xn h — Xi 

Vi 


)} 


dy 
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Let /3 > 0 be the constant in the statement of the lemma then for —h—M/m < Xi < 0, 
1 < i < n and < s < t/m, we have 


'Wn 


Qs{x,y)l(-h,h)-‘{y) dy 


> 




umn e dr. (56) 


Similarly, for 0 <Xi < h + M/m, 1 < i < n and s in the same range as above, we have 



Qs{x,y)l{-h,h)r^{y) dy 


> / Pi{y) TT HMY) e -V 2 M{tm)-y^)} dF. 

Juin) 


(57) 


Taking m large enough and noting that Pi(Y) is the probability density of a GUE matrix 
r, we see that both (56) and (57) can be made greater than (3 and hence completes the 
proof. □ 


Lemma 6.4. Let j3 be the constant in Lemma 6.3. Let t > 0, M >0 and h > 0 be 
such that {—h,h) C (—2M, 2M) and let Mn be the solution to (15) with initial data 
g = Then, there exists an mg := mQ{h, M,n,t) such that for all m > mg 


Mn{s,y) > ^l(-h-M/m,h+M/-m}^iy) for all t/2m < s< t/m and y e 


> 1—S{m), 


where S(jn) satisfies (1 — (5(m))™ —>■ 1 as m ^ oo. 


Proof. Let /3 be as in Lemma 6.3 and let M > 0, t > 0, h > 0 be given, then by Lemma 
6.3 there exist an mg = mg(/i, M, n, t) such that for all m > mg, all s £ [t/2m, t/m] and 
y G K" 

P —h—M/m,h+M/m)'^iy)' 

Since J„ is deterministic, we have 


Mn{s,y) < ^l(-h-M/m,h+M/m)’^iy) for some s G [t/2m,t/m] and y G R 
< p|^/„(s, y) < - ^l(-h-M/m,h+M/m)'^{y) for some s G [t/2m, t/m] and y G 


< 


sup \In{s,y)\ > - 

_yG{ —h — M/m,h-\-M/m)'^ 


lA 

:>|'Co 

1 

sup \Ln{s,y)\P 


\ ^ / 

s^[t/2m,t/m] 



_yG{ — h—M/m,h-\-M/m)^ 


<(iy\ 

sup |/„(s,y)|^’ 

V2/ 

{s,y)^[t/2m,t/m] X [—3M,3M]'^ 



( 58 ) 
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for all p > 2 by Chebychev’s inequality. We shall bound the final expectation. Fix 
a < 1/4 and 0 G (O, a — then since /„(0, y) = 0 for all y, we have 



Inis,y) 

P 

< E 


Inis,y) - Ini0,y) 

P 

sup 

{tjm)^ 


sup 

Q 


s£ [t/2m,t/m] 



s£ [t/2m,t/m] 



_yG[-3M,3M]” 




_yG[-3M,3M]” 




< E 


sup 

s,s' ^[0,t/m],s^s' 
yG[-3M,3M]^ 


In{s,y) - In{s',y) 
|s — s'l® 


p 


(59) 


Recall that Kolmogorov’s continuity criterion (see [RY99, Theorem 2.1]) states that for 
a stochastic process {X{t) : t G [0,T]'’*), if there exist strictly positive constants C, a 
and p with ap > d such that 

\\X{s)- X{t)\\p<C\s-t\°‘, foralls,tG [0,T]^, 


then X has a Holder continuous modification which satisfies for all 9 G [0, a — d/p), 


sup 

s^t 

s,tG[0,T]" 


|X(s)-X(t)| 


— t\^ 


< CT° 


2S+i2^/p 


1 _ 2 rf/p 2 -(“-®) ’ 


(60) 


Note that the right hand side of (60) is bounded for all p>2. 

By Proposition 5.3, for all p>2, there is a constant C := C{a,n) such that for all 
{s,y), {s',y') G [0,t/m] x [-3M,3M]", 

\\In{s,y)-In{s',y')\\p<Ccp sup ||M„(s, 2 /)||p(|s - s'|“ + I?/-(61) 

se[ 0 ,i/m] 

ye[-3M,3M]’^ 


Then by Kolmogorov’s continuity criterion, for p > (n + l)/a there is a constant K' := 
K'{a, M, TO, n, t) such that (59) is bounded by 

{K'rc; sup \\M„is,y)\\;<{4.K'^re^P"*/^, 

se[0,t/m] 

ye[-3M,3M]’^ 

for a constant A depending only on n, where to obtain the inequality we have used the 
moment bound (16) and the fact that 5 < 1, |erf(-)| < 1 and Cp < 2^. Furthermore, 
if TO > Too a t then tjm <\ and thus for such m we can, by the explicit bound on the 
right hand side (60), replace the constant K' in the previous display with a constant 
K := Kr(a, M, n). Consequently, for all p > (n + l)/a 


sup \In{s,y)\P 

sG [0,t/m] 

,yC[-3M,3M] 

< exp (- \-p\og{8 Ky/p) — p 01 og(m) 

\ TO 

Choose p = 4(n+ l)/a > (n+ l)/a and d = a/2 and for such p denote the exponenial in 
the last line above by 5{m), then for to large, 5(rn) ^ exp(— log(TO"+^)) and therefore 

(1 - <5(to))’” ^ (1 - ^ 1 , asm ^ 00 , 

for all n > 1 as required. □ 



< ('ll 
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We are now ready to prove the main result of this section. 


Proof of Theorem 6.2. By linearity is the solution to (15) with initial data 

gi — 92 and so it suffices to prove that y) > 0 for alH > 0 and y G R”] = 1, for 

g such that g > 0 and g{y) > 0 for some y G K" almost surely. 

We first consider the case when g is a continuous function such that g > 0 and 
g{y) > 0 for some y G M" so that one can find constants c > 0, d > 0 small enough 
such that g{x) > cHILi ^iyi-d,yi+d){x) for all x G K”. Without loss of generality, we 
can assume c = 1 and take y to be the origin for convenience. By the weak comparision 
principle (Lemma 6.1), we can therefore assume that the initial data is g{-) = l(-d,d)"(’)' 
From now on we drop the superscript g and just write M„(t, y). 

Let 'Y = /312 where /? is the constant in Lemma 6.3. Fix t > 0 and M > 0 such that 
(—d, d) C (—M, M). For fc = 1,..., m, define the events 

Ak := < Mn{s,y) > (d) for ad s G 



and for fc = 2,..., m the events 

Bi := |M„(t/2m,y) > 7 l(_d_M_d-HM)n(j/) for all y G M”| 

Bk ■= 

{k — l)t {2k — l)t 


Mn{s,y) > 7 l(-d Mk (y) for all s G 


m 


2m 


and y G 


We consider first the sets Ak. By Lemma 6.4, there is an mg such that for all m > toq 
there is a 5{m) such that 

IP[4li] > 1 — 6{m). 

Now assume that •-nAfe-i occurs. On the event we have M„((A:—l)t/m, y) > 
l^~^l(-d-M(k-i)/7n,d+M(k-i)/m)^{y) for all y G M" almost surely. Define a time shifted 
white noise by W^{s, y) = W{{k— l)t/m + s, y). Let M^{s, y) be the solution driven by 
the noise with initial data given by 7^“^l(-d-M(fe-i)/m,d-HM(fc-i)/m)"(d)- On the 
event Ak-i^ by the weak comparision principle, Mn{{k — l)t/m + s,y) > M^{s,y) for 
all s > 0 and y G R" almost surely. It is easy to see that M^{s,y) := y) 

is the solution to (15) with initial data l(-d-M(fc-i)/m,d-i-M(fe-i)/m)'*(d)- Lemma 6.4 
applied to M!f with h = d + M{k — l)/m shows that with the same mg and 5{-) as above 
that for all m > mg 


t t 
2m ’ m 


and y G 


P M„{s,y) > 7l(-d-i^,d-hi^)>*(y) for all s G 
and hence 

P M„"(s,y)>7H(_^_^ .d+^)^iy) for all s G 

By the above discussion, this implies that 

P[Afe|Ai n • • • n Ak-i\ > 1 — 5{m) for 1 < fc < TO. 


t t 
2 to ’ to 


and y G 


> 1 — i5(to). 


> 1 — 6{m). 


Now since Bi C Ai, P[i?i] > 1 — 6{m) and then proceeding in the same manner as 
before, we have 


¥[Bk\Bi n • • • n Bk-i] > 1 — S{m) for 1 < k < m. 
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Finally, by the union bound 


Q Afe n f] -Bfe 

,A:=1 


= 1 -: 


fl Afc u m -Bfc 


L \fc=i 

> 1- ( l-I 


n 

.fc=i 


\k^l 


- 1 -] 


.fe=l 


Since (1 — 5{m))'^ —>• 1 as m —>■ oo, we conclude that 
P[M„(s,y) > 0 for all s S (0,t] and y € [—M, M]"] > lim 


Pi n Pi Sfc 

.k—l k—1 


= 1 . 


Since t > 0 and M > 0 are arbitrary, this completes the proof in the case when the 
initial data g is a continuous function. 

We now prove the result for g satisfying the assumptions in Theorem 6.2(a). The 
idea is that after a small time r > 0, we are back in the situation above. We shall prove 
that for all T > 0, 


F[Mn{t, y) > 0 for alH > r and y S R"] = 1. (62) 

and since t is arbitrary this would imply the desired result. Let W'^{s, y) = W{t + s, y) 
be the time shifted white noise and let be the solution to (15) driven by the noise W'^ 
and with initial data Mn{T, •). The weak comparison principle shows that P[M„(t, y) > 
0 for alH > 0 and y G R"] = 1. We claim that P[M„(t, y) > 0 for some y] = 1 then 
since y i—>■ M„(r, y) is continuous, the strong comparison principle for continuous initial 
data proved above applied to the solution shows that F[M^{s,y) > 0 for all s > 
0 all y G R"] = 1 which proves (62). 

Therefore, it remains to prove the claim. Suppose the opposite is true, that is 
P[M„(t, y) = 0 for all y] > 0 and consider the solution M„(s, •) at time s < r. If 
Mn{s,y) > 0 for some y almost surely then the strong comparison principle for contin¬ 
uous initial data applies to show that M„(t, y) > 0 for all y almost surely. Hence, 
P[M„(s,y) = 0 for all y] > 0 for all 0 < s < r which implies that M„(0,-) = 0 
with strictly positive probability which is a contradiction. Thus, we must have that 
P[M„(t, y) = 0 for all y] = 0 which proves the claim. 

We now prove part (b) of the theorem; the everywhere strict positivity of M„(t, x, y). 
Fix T > 0 then the same argument as above together with Proposition 5.4 shows that 
P[M„(r, x,0) > 0 for all x] = 1. By the joint continuity of Mn, there exist random 
c = c(x) and d = d(x) strictly positive such that Mn{T,x,y) > cl(_£i,d)"(y) for all x, 
y G R” almost surely. For > 1 define the random set Bn := {x G R" : c(x) > 
l/N,d{x) > l/N}. Then Mn{T,x,y) > (l/Af)l(-i/Ar.i/Ar)"(2/) for aU y and all x G Bn 
almost surely. The strict positivity result proved above applied to the solution with 
initial data (l/Af)l(-i/Ar,i/A)’‘(2/) together with the weak comparision principle implies 
that 


P[i?Ar] := P[M„(t -I- s, X, y) > 0 for all s > 0 and y G R", x G Bn] = 1- 

By the joint continuity of M„, P[Ua=i = R"] = 1 and so P[nA=i = F[Mn{T + 
s, X, y) > 0 for all s > 0 and x, y G R"] = 1 as required. □ 
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